Supersolid Vortex Crystals in Rydberg-dressed Bose-Einstein Condensates 
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We study rotating quasi-two-dimensional Bose-Einstein-condensates, in which atoms are dressed 
to a highly excited Rydberg state. This leads to weak effective interactions that induce a transition 
to a mesoscopic supersolid state. Considering slow rotation, we determine its superfluidity using 
Quantum Monte-Carlo simulations as well as mean field calculations. For rapid rotation, the latter 
reveal an interesting competition between the supersolid crystal structure and the rotation-induced 
vortex lattice that gives rise to new phases, including arrays of mesoscopic vortex crystals. 
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Superfluidity, i.e., the frictionless flow of a liquid, is 
one of the most spectacular manifestations of quantum 
mechanical behavior on a macroscopic scale. Generally, 
superfluidity can be characterized from the response of 
a many-body system to a slow, externally imposed rota- 
tion. A classical fluid enclosed in a rotating vessel will 
be dragged along and eventually rotate with it. A super- 
fluid, on the other hand, will, due to a lack of viscosity, 
remain stationary, and, for sufficiently fast rotation, form 
quantized vortices that arrange on a regular lattice [U |2] . 
Superfluidity has long been speculated to occur even in 
solid states of matter, combining the seemingly antithet- 
ical qualities of crystalline order and non-dissipative flow 
[3]-[5 . Recent experimental evidence for such a peculiar 
supersolid phase in ^He crystals [6, 7 has sparked an in- 
tense debate about its physical origin, which, however, 
remains controversial [8 . 

On the other hand, ultracold gases have emerged as 
a powerful laboratory to study a diverse range of many- 
body problems, including discrete supersolid states in op- 
tical lattices [9 . A promising route to the observation of 
continuous supersolids has been recently laid out [lQl[TT], 
based on off-resonant dressing of atomic Bose-Einstein 
condensates (BECs) to high-lying Rydberg states [12]- 
[T6] . The effective atomic interactions resulting from such 
a Rydberg-dressing provide a clean realization of a simple 
model for supersolidity [17^ . Given the range of available 
techniques for accurate probing of BECs [18] and recent 
advances in the manipulation of cold Rydberg atoms [19 , 
this approach holds promise for the observation of super- 
solidity under well-defined and highly controllable condi- 
tions. 

In this Letter, we explore the possibility for creat- 
ing and detecting mesoscopic supersolids in confined, 
Rydberg-dressed BECs by probing their response to 
forced trap rotations. For infinitesimally small rota- 
tion frequencies, large-scale simulations demonstrate that 
crystalline order and superfluidity (see Figjl]) persist 
at well-accessible temperatures, enabling supersolid cre- 
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FIG. 1: (color online) (a) Schematics of the considered setup 
in which the groundstates (|^)) of condensed atoms are off- 
resonantly coupled to high- lying Rydberg states (|e)) with a 
large laser detuning A and a small Rabi frequency Q ^ A. 
The resulting effective interaction, shown in (b), can lead 
to the formation of mesoscopic supersolids in a quasi two- 
dimensional condensate, exemplarily shown in panels (c) and 
(d). The extended plateau of the histogram P{k) of /c-particle 
permutations, shown in (c), indicates a large superfluid frac- 
tion, while the density shown in (d) demonstrates crystal- 
lization. The results correspond to N — 3000 Rubidium 
atoms confined in a trap with C(;tr/27r = 125Hz and dressed 
to Rb(505'i/2) Rydberg states with realistic laser parameters 
of A/27r = 75MHz and VL/2ti = 2.2MHz. The depicted QMC 
results yield a large superfluid fraction /s = 0.51(4) for an 
experimentally accessible temperature of T = 42nK. 



at ion in cold atom experiments. Combining the results 
from first-principle Quantum Monte Carlo simulations 
and mean field calculations, we reveal a universal be- 
havior of the crystalline ordering and the superfluidity 
that enables simple estimates of the system behavior over 
a wide range of experimentally relevant conditions. At 
higher frequencies, the rotation induces drastic structural 
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changes due to an interesting competition between the 
vortex lattice, controlled by short-range collisions, and 
the supersolid crystal arising from the finite-range Ryd- 
berg interactions. 

A Rydberg-dressed BEC is a gas of Bose-condensed 
atoms whose ground state \g) is far off-resonantly cou- 
pled to a highly excited Rydberg state |e) with a Rabi 
frequency 1] <C A much less than the corresponding laser 
detuning A fsee Fie; jlja)] [T0HI2. Specifically, we con- 
sider two-photon coupling to nS'1/2 Rydberg states of al- 
kalines which leads to an effective interaction [10] 



W{r) 



(1) 



between the dressed ground state atoms, where Cq = 
{^YCq, Rc = (C6/2nA)i/6 and Ce is the strength of 
the van-der-Waals interaction between Rydberg atoms 
[20] . At large distances r ^ R^^ the interaction resem- 
bles the Rydberg- Rydberg atom interaction Ce/r^, sup- 
pressed by a factor since only a small Rydberg 
state fraction ((7/2 A) ^ is admixed to the ground state. 
At small distances the dipole blockade [21 prevents two- 
atom dressing, such that W{r) saturates to a value of 
(^) hTt^ independently of the addressed Rydberg state 
[see Fig 1(b)]. 

We consider quasi two-dimensional condensates [22, 
[23] . that are strongly confined along the z-direction. 
Transversely, the atoms are weakly confined in a har- 
monic trap with a frequency c^tr- Scaling lengths by the 

^ (m is the atomic mass) and 



oscillator length / 



frequencies by cjtr, we obtain the following Hamiltonian 
in the rotating frame 



i<3 



+ 7<5(fi-f,), (2) 



where 7 = y^na^/lz is the effective strength of the con- 
tact interaction due to s-wave collisions with scattering 
length tts, a = mCQ/ ih^l^)^ Vc = Rc/h ^ is the rotation 
frequency in units of the trapping frequency, = Vi/l is 
the scaled position of the ith atom and Lz^i denotes the 
corresponding angular momentum operator. 

Focussing on the effects of the Rydberg interactions we 
first assume 7 = and investigate the stationary {Vt = 0) 
equilibrium thermodynamics of H by Quantum Monte 
Carlo (QMC) simulations [2j based on the continuous- 
space Worm algorithm [25]. Figjljd) shows the com- 
puted density profile for a large number {N = 3000) 
of Rubidium atoms at a temperature T = 42nK cou- 
pled to Rb(505'i/2) Rydberg states with Q/27r = 2.2MHz 
and A/27r = 75MHz. The formation of a self-assembled 
droplet crystal is evident. Note that the experimentally 



relevant parameters correspond to very weak interactions 
k-^W{r) < k-^hn^/{2A)^ - 10-^nK <C T. In this weak 
coupling limit, crystallization, nevertheless, occurs due 
to the presence of many atoms in each blockade sphere, 
forming a solid of liquid droplets with collectively en- 
hanced interactions. On the other hand, the average 
number of Rydberg excitations in each droplet is only 
~ 0.04, which diminishes Rydberg state decay [TOl [13] 
and renders many-body effects, that tend to counteract 
the binary interaction ([T]) [26 , unimportant. 

In addition, the QMC simulation yields a large su- 
perfiuid fraction of /s ~ 0.5, showing that the observed 
mesoscopic crystal is indeed a supersolid. /s corresponds 
to the fraction of atoms that decouples from a vanishingly 
small trap rotation, such that = 1 — Iqm/Ici where Ic 
and Iqm are the classical and the actual (i.e. quantum 
mechanical) moments of inertia, that can be computed 
within QMC [27]. The finite value of fs is also indicated 
in Figjljc) by the extended plateau of the probability 
P{k) that k particles exchange within the system. P{k) 
is finite for all /c < A/" in a superfiuid, but drops to zero 
in a classical crystal [28 . 

The large superfiuid fraction arises from the extended 
inner plateau of the interaction potential ([T]) which can 
accommodate a large number of atoms that maintain a 
sizable superfiuid fraction, even in the defect-free droplet 
crystals discussed in this work. As a result, the superfiu- 
idity does not critically depend on the size and specific 
structure of the mesoscopic crystals. This behavior is in 
contrast to mesoscopic assemblies with pure power-law 
potentials [29 where the superfiuidity depends strongly 
on the cluster geometry and vanishes in the bulk-limit of 
infinite defect-free crystals [30] . 

The temperature dependence of /s is shown in Figj2|a) 
for different atom numbers N and interaction strengths 
a, but with the product aN = const, held constant. Ap- 
parently, superfiuidity extends to higher temperatures for 
larger values of which is readily understood from sim- 
ple arguments. As will be explained below, the crystal 
structure remains unchanged for aN = const., such that 
increasing N merely increases the number of atoms in 
each droplet. As a result, the system becomes less sus- 
ceptible to fiuctuations, which facilitates exchange be- 
tween the mesoscopic crystal sites and, hence, leads to 
larger /g. Upon rescaling the temperature with A/", all 
data collapse on a single curve [inset of Figj2|a)], which 
is consistent with the linear density scaling of the crit- 
ical temperature for the Berezinskii-Kosterlitz-Thouless 
transition ^31] in the weak-coupling limit [32 • This uni- 
versal behavior of the crystalline order and superfiuidity 
with respect to aN and T/N can be used to predict the 
system behavior for a wider range of parameters. For 
instance, at a typical interaction range Rc = 2/im and 
a large number of A" = 10^ Rb atoms, the calculations 
of Figj2] imply a crystal structure as in Figj2|e) with a 
superfiuid fraction of /g ~ 0.4, both persisting at a large 
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FIG. 2: (color online) (a) Temperature dependence of the 
superfluid fraction for different particle numbers and aN = 
24500, 7 = and Tc = 2.65. The horizontal arrow indicates 
the T = MF prediction [eq.(|3|]. As shown in the inset, 
all data points collapse on a single curve upon scaling T by 
N. The vertical arrows indicate the temperatures used in (d) 
and (e) . (b) Comparison of the radial density p and superfluid 
density ps for N = 400. Panels (c)-(e) show the corresponding 
ground state densities obtained from (a) the MF approxima- 
tion eqs.Q and (|3| and (b,c) the QMC simulations. 



temperature of T ^ 350nK. 

The local superfluid density ps is also closely linked to 
slow trap rotations, and can be obtained within QMC 
from a local area-estimator [33 of the non-classical ro- 
tational inertia. In his seminal paper on supersolidity 
in Helium, Leggett [5 proposed a far simpler approach 
to estimate the superfluid density by considering the re- 
sponse of a liquid in a slowly rotating ring-shaped con- 
tainer. Within a mean field (MF) approximation, sim- 
ple phase matching arguments yield the radial superfluid 
density [5j 



Ps{r) 



27r 



Jdipp{r)- 



(3) 



solely in terms of the single-atom density p(r), integrated 
over the polar angle cp. This expression provides an up- 
per bound of the exact superfluid density Ps{t) and offers 
an intuitive picture for the superfluid properties of the 
found droplet crystals, suggesting that the maintenance 
of a superfluid flow is foremost hampered by regions of 
low density O |34l |35] . With increasing interactions the 
density modulations become more pronounced such that 
Ps{t) eventually vanishes, since the atoms can not pass 
through extended regions where p(r, (p) ^ 0. In this 
strong-coupling limit, that goes beyond the MF approx- 
imation, QMC simulations show that the system enters 
a globally insulating phase [11]. Likewise, eq.([3| gives 
/)s(r) = p{r) in the non- modulated phase, i.e. a perfect 
superfluid. 



In order to test this intuitive picture we also performed 
MF calculations, based on the corresponding nonlocal 
Gross-Pitaevskii equation (GPE) 



idt'ip{Y) 



(4) 



-fN\tlj{r)\^ ^ aN 



+ r 



where tp{r) denotes the normalized (J |?/^pdf = 1) con- 
densate wavefunction and yields the single atom density 
p{r) = N\tjj{r)\'^. Eq.Q predicts that within MF the 
crystal structure depends only on the product aN^ as 
found above from the QMC simulations. At very low 
temperatures, the MF calculations yield virtually the 
same crystal structure as the QMC results [cf. Figj2|c) 
and (d)]. In this zero-temperature limit the superfluid- 
ity approaches a constant value independent of the atom 
number N [see Figj2|a)], with aN held constant. This 
behavior can also be understood from eqs.(|3| and Q, 
which depend solely on the product aN. Indeed, the 
radial atomic density and superfluid densities obtained 
from eqs.Q and Q closely reproduce the QMC results, 
demonstrating that eq.(|3| not only provides an upper 
bound O [35] of the superfluidity, but yields an accurate 
description of the structural and superfluid properties for 
the present realization of a supersolid. 

Having established the predictive power of the MF ap- 
proximation for T ^ 0, we can apply eq.Q to rotating 
assemblies. Previous MF studies of rotating dipolar con- 
densates [36] have revealed a range of different vortex 
lattices as a function of the dipolar interaction strength. 
In the present case, the ground state depends on four 
parameters Tc, jN and aN^ which makes a complete 
parameter scan prohibitively demanding. Being primar- 
ily interested in the effects of the Rydberg interactions 
and focussing on the large- regime, we set = 10"^ 
[37j . In order to scan the remaining parameter space, we 
use the following procedure: For a given set of Q and Vc 
we first calculate the ground state for a = via complex 
time evolution, which agrees with the results [38 for con- 
tact interactions. Subsequently, we slowly increase aN 
and follow the complex time evolution [TOl |38] of the 
BEC. This yields the "phase diagrams" depicted in Fig[3| 

In the region of small aN and Q (SF in Figj3| we 
find a simple superfluid, with p^ = p and an unstruc- 
tured density profile. Upon increasing the rotation fre- 
quency above a critical value, vortices start to form and 
eventually arrange on a regular lattice, whose spacing 
decreases for increasing rotation frequency and interac- 
tion strengths (VL in Figj3|. For small values of aN 
the nonlocal character of the corresponding interaction 
in eq.Q is of minor importance, such that the total in- 
teraction can be described via an effective contact term 
7^ = 7+27r^a/ (3^/^r^). In this regime, we indeed find or- 
dinary vortex lattices determined by ^y'N [see Figsjsj^a)^, 
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FIG. 3: (color online) "Phase diagrams" for = 10^ and (a) rc — 1, (b) Tc — 3, (c) Tc — 5. The insets display ground 
state density profiles at the indicated parameters, showing only the inner trap region for better visibility, where dark shading 
indicates high density and vice versa. The dimension of each inset is indicated by the vertical bar corresponding to a length 
of 5/. The horizontal arrows at Q = show the supersolid transition obtained from LDA estimates of the roton instability as 
described in the text. 



(b) i and (c)i]. With increasing aN the finite-range of 
the dressing-induced interaction becomes important and 
changes the symmetry of the vortex lattice. For small 
[Figjsja)] the system is already close to the bulk limit, 
where one finds a transition to a honeycomb density pat- 
terns. For large finite size effects become dominant 
and lead to the formation of concentric rings superim- 
posed on the underlying vortex lattice [Figj3];^c)n]. 

The finite system size also affects the crossover to su- 
persolid states [39]. In the bulk limit, the critical inter- 
action strength 0^20^ can be linked to a roton instability 
of the unmodulated groundstate p^, and, therefore, de- 
creases with increasing rotation frequency (at larger ^ 
it increases again due to the decreasing density caused 
by the centrifugal potential). For small [Figjsja)], the 
critical value at 1] = can, thus, be estimated from a 
local density approximation (LDA) of the roton instabil- 
ity, which gives r~^aN\'il){0)\^ = 31.9 + 6.I7 and for the 
parameters of Figjsja), r~^aN = 9.3 x 10^. For larger 
Vc the supersolid transition is preceded by the aforemen- 
tioned ring formation. Assuming that the transition is 
initiated by a one-dimensional roton instability of the 
innermost ring one obtains r~^ai-DN = 6.7 x 10^ and 
r~^ai-[:,N = 7.2 x 10^ for the parameters of Figjsjb) and 

(c) , respectively, in good agreement with the numerics. 
The supersolid phase (SS in Figj3| retains the vortices. 

They are, however, pushed into the low-density regions 
and form a honeycomb vortex lattice that does not af- 
fect the structure of the supersolid droplet crystal (see 
Figsj3jb)in and (c)m). At larger Tc, on the other hand, 
there is an interesting competition of length scales be- 
tween the supersolid crystal and the vortex lattice. For 
sufficiently high rotation frequencies the vortex density, 
and, hence, the vortex-vortex interaction, exceeds a crit- 



ical value at which it becomes energetically favorable 
to form vortices inside the superfluid droplets [SSV in 
Figjsjc)]. As shown in Fig jsj^c) ii;, the crystal structure 
of this additional vortex lattice is imposed by the trian- 
gular supersolid. The minimum rotation frequency, to 
support these states increases as Tc is reduced. Since 
the confinement ceases for Q > 1, this implies a mini- 
mum interaction range r^^^{jN) in order to observe the 
SSV states for a given contact interaction 7A^. For larger 
Tc the on-site number of vortices successively increases, 
which eventually form a mesoscopic crystal of small vor- 
tex crystals as shown in Figj4j 

In this work, we demonstrated novel vortex transitions 
in a realistic model of supersolids in the regime of weak 
interactions and large atom numbers per blockade sphere, 
which is most directly accessible to experiments. Extend- 
ing the present study to the regime of strong coupling, 
future work may address how strongly correlated super- 
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FIG. 4: Density for Q = 0.8 and (a) = 7, aN/rt = 
7.5 X 10^ and (b) rc = 9, aN/rt = 9 x 10^. 
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solids can be stabilized in finite-temperature cold atom 
experiments, which will more closely resemble the phys- 
ical mechanism thought to underlie super solidity in He. 
Along these lines, the realization of the found vortex crys- 
tals in the regime of strong interactions and small particle 
numbers will open the way towards studying quantum 
Hall-like phenomena in rapidly rotating self-assembled 
structures. This will require a beyond-MF description 
of the correlated many-body dynamics induced by rapid 
rotation. 
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